Circuit QED with a Flux Qubit Strongly Coupled 
to a Coplanar Transmission Line Resonator 
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We propose a scheme for circuit quantum electrodynamics with a superconducting flux-qubit 
coupled to a high-Q coplanar resonator. Assuming realistic circuit parameters we predict that it 
is possible to reach the strong coupling regime. Routes to metrological applications, such as single 
photon generation and quantum non-demolition measurements are discussed. 



I. INTRODUCTION 

Until a few years ago it was an open question whether 
true quantum effects such as quantum entanglement 
would ever be observed in a man-made macroscopic elec- 
tronic device. However, over the past decade, quantum 
coherence has been demonstrated in a variety of macro- 
scopic systems, including superconducting circuitsii^i^i^. 
Many of these experiments drew inspiration from the pi- 
oneering work on atomic qubits that took place a decade 
earlier—. As the fields of atomic physics and quantum 
optics continue to advance, it makes sense to continue to 
look to them for guidance. 

The universal nature of quantum mechanics is greatly 
to our advantage, in that the terminology and method- 
ology apply as well to macroscopic as to microscopic 
systems. This allows well-known results from atomic 
physics and quantum optics to be used to plan and pre- 
dict the outcome of experiments on solid state devices. 
Recently, such techniques have been applied with great 
success to implement a number of ideas such as quan- 
tum state tomography^, Mach-Zehndcr interferometry^ 
and sideband cooling*. 

This approach has also been very successful in the field 
of cavity quantum electrodynamics (CQED) [see Ref.— 
for an introduction]. In CQED an atomic 2-level system 
(i.e. a qubit) is made to interact with a high-finesse op- 
tical cavity with a coupling energy Tig. Provided that 
the relaxation rates 7 of the qubit and k of the cavity 
field are smaller than g (known as the strong coupling 
criterion), it is possible to observe a coherent exchange 
of energy between the qubit and the cavity field. The 
resulting entangled states can be detected spectroscopi- 
cally. Recently, Schoelkopf and co-workersi2iii achieved 
strong coupling in a macroscopic circuit comprising a su- 
perconducting charge qubit and a coplanar transmission 
line resonator. This new field is known as circuit-QED, 
and has many potential applications such as the gener- 
ation and detection of single microwave photons. More 
recently strong coupling was observed in experiments on 
photonic crystalsiS and quantum dots^J,. 

Much of the work to date on qubit-cavity sys- 
tems has been focused on methods for making quan- 
tum non-demolition (QND) measurements of the qubit 
state. QND schemes have been used to read-out single 



qubitsiSii^ and also to measure the photon number in 
the cavity^^. Several ways of reading out qubits using 
low-Q cavities/resonators have also been reported^®' 
but these do not allow the strong coupling regime to be 
reached. 

The benefits of superconducting circuit-QED systems 
over atomic systems are twofold. Firstly, the qubit en- 
ergy can be tuned by varying the external magnetic field, 
enabling control over the qubit-cavity interaction. Sec- 
ondly, qubit parameters such as the level separation at 
the degeneracy point can be engineered through appro- 
priate circuit design. The main advantage of flux qubits 
over other types of superconducting qubits is that they 
are less susceptible to fluctuations in the background 
charge and the associated noise. This makes them less 
prone to decoherence, and therefore easier to manipulate 
in a deterministic way. 

In section II we summarise the established background 
theory in the context of our proposed system; in III we 
show that it is possible to achieve strong coupling be- 
tween a superconducting flux qubit^i^ and a high-Q copla- 
nar transmission line resonator; in IV we present com- 
puter simulations of the resonator response; and in V we 
propose a scheme for producing single photons on de- 
mand. 



II. THEORY OF CIRCUIT QED WITH A FLUX 
QUBIT 

In this section we analyze the processes which occur 
when a superconducting qubit is coupled to a supercon- 
ducting coplanar transmission line resonator, as shown in 
Fig. [TJ The effect of coupling a fiux qubit to a resonator 
has previously been experimentally demonstratedi^ but 
the quality factor of the resonator was too small to fulfill 
the strong coupling criteria. 

Two types of flux qubits will be discussed - an RF 
SQUID, which consists of a single Josephson junction in 
a superconducting loop, and a persistent current qubit 
(PCQ), which has three junctions in the loop, one of 
which is smaller than the others by a factor a^^ . Both 
need to be biased by an external magnetic flux <^x , which 
tunes the energy level separation through an anticross- 
ing at = 0.5$o- The resonator is a coplanar trans- 
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FIG. 1: a) Sketch of a typical coplanar waveguide resonator 
of length l=A/2 « 11 mm. Shown is also how the qubit can be 
placed in between the centre conductor and the ground plane 
of the waveguide, b) Schematic diagram of a superconducting 
qubit coupled to a coplanar transmission line resonator. (A) 
Persistent current qubit. (B) RF SQUID. M is the mutual 
inductance between the qubit and resonator, and $x is the 
magnetic flux threading the qubit loop. 

mission line with inductance L and capacitance C, which 
is weakly coupled to external transmission lines via cou- 
pHng capacitors Cc- The Hamiltonian H of the complete 
system is 

H = Hg + Hr + Hg + Hi + HE. (1) 

Hq describes the qubit, Hr the resonator, and Hg the 
interaction between them. Hi denotes the interaction 
of the system with a periodic drive field. Finally, He 
describes the interaction of the resonator with its envi- 
ronment, resulting in the loss of photons to the external 
transmission lines and the interaction of the qubit with 
its environment, resulting in spontaneous decay from the 
excited state to the ground state. 

The qubit Hamiltonian Hg is given by the expression 
(?i/2)(— cCTz — Sax) where cr^ and are Pauli spin matri- 
ces, 6 is the the level repulsion, e = {2Ip/h){^x ~ "J'o/S), 
and Ip is the persistent current. In the case of an RF 
SQUID suitable for operation as a qubit, Ip is roughly 
equal to half the critical current of the single junction, 
whereas, for the persistent current qubit, Ip is approxi- 
mately equal to the critical current of the smallest of the 
three junctions. If the qubit is operated at or near the 
degeneracy point, Hq can be expressed more simply by 
transforming to the basis in which the ground state | |) 



and excited state | t) correspond to symmetric and anti- 
symmetric superpositions of clockwise and anti-clockwise 
persistent currents. This yields 

^, = ^(IT>(T|-U>ai) = ^-., (2) 

where the level separation is Tuoq, luq being the qubit 
Larmor frequency Ve^ + . 

Assuming the resonator supports only a single mode 
of the electromagnetic field, its Hamiltonian is given by 

Hr — hujr ^a^a -I- , (3) 

where a^{a) is the creation(annihilation) operator which 
creates (destroys) a single photon in the cavity. The 
eigenstates of the resonator described by this Hamilto- 
nian are Fock states |0) . . . |n) . . . with n photons. The 
single mode condition corresponds to a harmonic oscilla- 
tor with the energy levels at TiLJr{n -1-1/2) and the zero- 
point energy hiOrl^- 

The interaction between the radiation and the qubit 
is described as the dipole interaction Hg = —jl ■ B be- 
tween a magnetic moment /i of the persistent current 
circulating in the qubit loop and the magnetic field B in 
the resonator. Introducing quantization, this term in the 
Hamiltonian can be written in the form 

Hg = Tig{a)a^ + cr^a), (4) 

where a~^{a~) is the raising (lowering) operator for the 
qubit. This expression is valid in the so-called non- 
dispersive regime where ujr « wo^. The constant g char- 
acterizes the qubit-photon interaction strength and the 
expression in the brackets describes the process whereby 
the qubit can be excited by absorbing a photon, or a 
photon can be generated at the expense of de-exciting 
the qubit into its ground state. In the next section we 
shall explicitly calculate the dipole coupling strength g 
for specific designs of the qubit and the cavity. 

The interaction of the coupled system with an external 
classical drive field can be seen as a periodic exchange of 
photons between the resonator and the driving field: 

i?/ = e(e"'"*a^ + e"^*a), (5) 

where ^ is the drive amplitude. One can also drive the 
qubit directly using a separate control line leading to 
terms ^'(e~"^*o'+ + e*"*(T^), but we shall not consider 
this case any further. 

Adding the terms (2)-(5) together we arrive at the 
driven Jaynes-Cummings Hamiltonian^^ 

H — ~^'^z + iT-LOr ^a^a + + hg{a) + a) (6) 

-f ^(e-*'^*at -He*'^*a). 

This Hamiltonian can be used to write down a mas- 
ter equation which completely describes the dynamics of 
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the system. All the parameters which define the system 
can be conveniently written in units of angular frequency 
- a convention which we will follow here. If the cav- 
ity is weakly coupled to an already weak drive field one 
can achieve a regime when only two lower Fock states of 
the resonator are relevant. Within the picture described 
above we make one two- level system, the qubit, interact 
with another two-level system, the resonator. When the 
qubit is detuned from the resonator the eigenstates of 
the coupled system can be written: |0 |), |0 t)) |1 i) 
and |1 t)j where the number represents a Fock state of 
the resonator and the arrow represents the qubit state. 
However, when the qubit is brought into resonance, Hg 
couples the states |0 t) and |1 |) and lifts their degen- 
eracy. The system will oscillate between the states |0 t) 
and |1 I) at a frequency 17 = 25 , known as the vac- 
uum Rabi frequency^!, giving rise to a splitting^^ of the 
central peak as shown in Fig. [21 One can visualize this 
as a cycle in which the resonator and qubit continuously 
exchange an amount of energy equal to one photon. As 
the drive amplitude increases it will start to perturb the 
system. This leads to a set of states that are shifted by 
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Hence, the effect of the drive is to reduce the Rabi fre- 
quency. 

The coupling strength g can be determined experimen- 
tally by making a spectroscopic measurement of the split- 
ting AEi. The drive amplitude should be reduced until 
the splitting reaches its maximum value, where it is equal 
to 2gh. 
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FIG. 2: Energy levels of the coupled qubit-cavity system. On 
the left of the diagram the qubit is far detuned from the cavity. 
As we move from left to right, the magnetic flux threading the 
qubit loop is increased, tuning the qubit transition frequency 
into resonance with the cavity. On the right of the diagram, 
the qubit is once again far detuned. 



The above discussion covers the resonant regime, 
where the qubit is tuned into resonance with the cav- 
ity. In contrast, when the detuning A = loq — Wr is large, 
such that g/A <C 1, a dispersive Stark shift pulls the cav- 
ity frequency by ±g^/A. This so-called dispersive regime 
can be used to perform quantum non-demolition mea- 
surements of the qubit stated. 

The effects of the environment on the system are taken 
into account by the term He- There are three types of 
damping that need to be considered: 

• Photons leak out of the cavity at a rate k = uJr/Q, 
where Q is the cavity quality factor. 



• The qubit relaxes at a rate 7 
the energy relaxation time. 



1/Ti, where Ti is 



• Pure dephasing of the qubit at a rate 7^ = 1 /T^ = 
I/T2 — 2/Ti, where T2 is the dephasing time. 

Dephasing plays a larger role in solid state systems than 
atomic systems, due to the stronger interaction of solid 
state qubits with their environments. In the absence of 
pure dephasing we would have T2 = 2Ti but in real sys- 
tems T2 is frequently much shorter than that, indicating 
the need to take pure dephasing into account. 



III. STRONG COUPLING WITH A FLUX 
QUBIT 

Below we estimate the coupling strength g using a 
semi-classical approach. We treat the flux qubit as 
a magnetic dipole and assume that it is placed at a 
magnetic field antinode of the resonator. The coupling 
strength is given by 5 = /iBorms/^, where Borms is the 
zero-point root mean square magnetic field generated by 
the current fluctuations at the antinode of the resonator. 
The magnetic dipole moment of the qubit is given by 
fi — IpA, where Ip is the persistent current flowing 
around the loop and A is the loop area. We can esti- 
mate Borms by considering the zero point energy of the 
resonator, fujJr/2. This energy cycles continuously be- 
tween inductive and capacitive components. The mag- 
netic field is determined by the inductive component, 
LI{tY /2, where L is the total equivalent inductance of 
the resonator near resonance and I{t) is the instanta- 
neous current. At the moment when the energy is purely 
inductive, we have (l/2)L/^j^j. — {l/2)hu!r- Since I(t) 
undergoes sinusoidal oscillation we therefore have that 
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Assuming that current flows in thin strips (whose width 
is determined by the superconducting penetration depth) 
at the edges of the centre conductor and ground plane, 
the field at the antinode of the fundamental mode is ap- 
proximately given by 
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where /xq is the permeabiHty of free space and r is half the 
width of the gap between the centre conductor and the 
ground plane (we assume the qubit is placed at the centre 
of the gap). Therefore, the coupling strength between the 
qubit and resonator is given approximately by 

By inserting realistic values for the parameters in the 
above equation we can obtain an estimate of g. 

First, we choose the fundamental frequency of the res- 
onator. It is convenient to choose a value that lies within 
the range 4-8 GHz, as this is well within the design scope 
of both the qubit and resonator, and can be accessed 
with commercial microwave sources and components. We 
choose LOrjl'K = 6 GHz. With a centre conductor of width 
~ 10 /im and a gap of width ~ 5 /im, it is possible to 
achieve a total inductance L ^ 2 nH for a resonator op- 
erated at its fundamental frequency.. 

Next, we choose Ip such that the transition frequency 
tJo of the qubit at the degeneracy point is slightly less 
than that of the resonator. This will enable us to tune 
the qubit in and out of resonance with the resonator by 
changing the external flux threading the qubit loop. 
For the 3-junction persistent current qubit having two 
junctions of critical current Ic = 800 nA and junction 
capacitance C = 4 fF, and one junction of critical cur- 
rent a/c, where a = 0.72, we get loqII-k — 4.9 GHz at 
the degeneracy point and /p k, 580 nA. These parameters 
were obtained by solving the Schroedinger equation nu- 
merically. The transition frequency cjq does not depend 
on the area A of the qubit loop, provided that the loop 
inductance remains small compared with the Josephson 
inductance. However, we note that the larger the loop 
area, the greater the (undesired) coupling to the environ- 
ment. Here we choose a value A w 8 [lvc? . 

With the above parameters we obtain 5/27r « 35 MHz. 
We now compare this with the rate of photon loss from 
the resonator k and the relaxation rate of the qubit 7. 
When g > k, 7, the coupled system is able to undergo 
many cycles (~ 2g / (k+j)) of vacuum Rabi oscillation be- 
fore losing coherence. This is important for applications 
such as single microwave photon generation. The photon 
loss rate from the resonator is given by /c/27r = ujr/Q, 
where Q is the loaded quality factor of the resonator. It 
is possible to design a resonator with Q = 10^, yield- 
ing K/27r « 0.1 MHz. The relaxation rate of the qubit 
is given by 7 = 27: /Ti. Taking Ti ~ 1 ^s2., we obtain 
7/27r « 1 MHz. Naturally, the values for Ti and T2 for 
a real system can not be predicted with any accuracy 
and will depend on the experimental conditions. How- 
ever, based on the data available in the literatureS^i^ we 
believe that the aforementioned values are reasonable. 
Hence, it is clear that our estimated value of g for the 
persistent current qubit should satisfy the strong cou- 
pling criterion. 

For an RF SQUID with critical current Ic — 10 ^P^, 
area A = 64 /im^, loop inductance L squid = 35 pH and 



junction capacitance C = 50 fF, we obtain a transition 
frequency 10^121: = 4.6 GHz. In contrast to the persis- 
tent current qubit, the area of the RF SQUID does affect 
the transition frequency, via the loop inductance. It is 
difficult to reduce the area further than the value we have 
chosen, as this necessitates increasing the critical current 
and decreasing the junction capacitance, which becomes 
increasingly difhcult to achieve in practice. Close to the 
degeneracy point, the persistent current Ip in the above 
SQUID is expected to be about 5 /xA. Combined with the 
increased loop area, this is likely to lead to an even larger 
coupling strength g than predicted for the persistent cur- 
rent qubit (unless the SQUID is displaced significantly 
from the antinode of the resonator) but at the same time 
make qubit more susceptible to noise. The fact that the 
loop area of the 3-junction persistent current qubit can 
be made small enough to render it relatively insensitive 
to flux noise is one reason why it has been so successfully 
used by e.g Mooij and co-workers^ii^. 



IV. NUMERICAL SIMULATION OF THE 
SPECTRUM UNDER MICROWAVE 
EXCITATION 

Having shown that it is possible to reach the strong 
coupling regime with a flux qubit coupled to a coplanar 
resonator, we now simulate the results of a spectroscopic 
experiment to measure g. The experiment would involve 
driving the coupled system with an external microwave 
field whose frequency uji would be swept through the res- 
onance of the coupled system. The most straightforward 
way to probe the response of the system would be to 
use what effectively amounts to a standard microwave 
transmission ( S12 ) measurement of the cavity. The ex- 
periment would be done by starting with the qubit far 
detuned from the resonator, then stepping the external 
magnetic flux to tune the qubit through the cavity reso- 
nance. This type of experiment would allow us to record 
the output power of the resonator as a function of qubit 
Larmor frequency. 

These simulations were performed by solving the mas- 
ter equation using a Liouvillian with the Hamiltonian ([7]) 
and three collapse (Lindblad) operators which account 
for the decay and dephasing of the qubit and the cavity 
at the aforementioned rates k,7 and 7^ (i.e. the effects of 
He)- a brief description of the formalism can be found 
in appendix [BJ All simulations were performed using the 
"Quantum Optics Toolbox" developed by Tan3^. Note 
that all the figures show the spectrum of the intracavity 
field (see appendix [Cl for the definition and a description 
of how it is calculated). 

In the regime where the qubit Larmor frequency is 
very far detuned from the cavity (i.e. when even disper- 
sive effects are negligible) the qubit and the cavity are 
effectively decoupled and no exchange of energy can take 
place. If the system is probed by measuring the response 
of the cavity, a single peak located at the bare resonance 
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FIG. 3: The spectrum far from resonance (where the only 
effect of the qubit is to broaden the resonance) and at zero 
de-tuning for small drive (mean steady-state photon number 
(n) < 10"'^), the latter giving rise to a Rabi splitting 2g. 
Shown is also the spectrum at large ((n) ~ 8) drive amplitudes 
{dashed line). 



frequency ujr will be seen. 

If instead the qubit is tuned exactly on resonance (A — 
0) the effects of the coupling become clearly visible. Now, 
there are two peaks located at i ^^/2 = i 5 as can 
be seen in fig. [3l In between these two extremes there 
is a gradual change from a single- to a double-peaked 
spectrum where the splitting is approximately g^/A as 
can be seen in the left picture of fig. U) In this case the 
pure dephasing rate 7^ was set to zero. The result is a 
"diamond-shaped" picture with a maximum splitting of 
2g at zero detuning. 

Far from resonance we can identify the four branches as 
being associated with the states |0 t) and |1 |) above and 
below the resonance. Exactly on resonance the system 
is in a superposition of states |0 t) i |1 i)- This is in 
agreement with the diagram shown in fig. [2] Note that 
e.g. an interferometric measurement method^- must be 
used in order to be able to observe the whole spectrum, 
a simpler transmission experiments would only see one 
peak off-resonance since such a measurement only records 
transitions between photon states where the final state 
can emit a photon; in this case |0 t) ^ |1 i)- 

The right picture in fig. [4] shows the spectrum when 
pure dephasing is introduced to the model. The result is 
somewhat more complicated in that the spectral weights 
are asymmetric with respect to Wr — wi — when the 
qubit is detuned from resonance. The reason for this 
is that dephasing leads to a loss of coherence, meaning 
that the qubit tends to stay in its ground state and the 
coupling between the states |0 t) and |1 |) is reduced, 
the system is therefore no longer in a superposition of 
those states. Whence, states with zero photons in the 
cavity are effectively "decoupled" from the onc-photon 
states. Since only states which allow for (at least) one 
photon in the cavity can be measured it follows that only 
the two branches (one above and one below resonance) 
that (approximately) correspond to ±|1 |) will be clearly 
visible. Note, however, that despite the loss of coherence 
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FIG. 4: (Color online) Power spectrum of the coupled qubit- 
resonator system as a function of qubit detuning in the strong 
coupling limit. The frequency of the drive field uji is held at 
the resonance frequency uJr of the bare resonator, while the 
Larmor frequency u>o of the qubit is tuned by changing the 
external magnetic field. Left: Spectrum in the absence of 
pure dephasing. Right: Adding a pure dephasing channel to 
the dissipation results in an asymmetric spectrum, here the 
pure dephasing rate 7^ is 9.5 MHz. Shown is also one branch 
of the expression ig^/A (white dashed line). The following 
parameters were used in the simulations (in GHz) : lj,. — uji = 
6, 5 = 0.035, K = 0.004, 7 = 0.001 and ^/R = 0.25k. 

an on-resonance measurement would still show two peaks 
separated by 2g, i.e. exactly on resonance this situation is 
effectively indistinguishable from the more coherent case 
even when the full spectrum is measured. 

While there are no bound states in the limit of strong 
driving S, > g/2 a continuum of states still exists giv- 
ing rise to complex spectra (dashed line in fig. The 
structure is reminiscent to the so-called "Mollow" peaks, 
well known in atomic physics from e.g. fluorescence spec- 
troscopy. However, in the latter case the peaks are the 
result of strong driving of the atom (qubit) whereas in 
this simulation the cavity is being driven so the similarity 
is somewhat superficial. When as in this case the cavity 
is driven so strongly that there are on average several 
photons in the cavity, we see both a drive induced shift 
of the position of the sidebands and a reappearance of 
the central peak. Note that since it is difficult to directly 
relate the parameter ^ to the power output from a mi- 
crowave generator, care must be taken not to drive the 
system inadvertently into this regime. 

For the persistent current qubit considered here, a de- 
tuning of ±0.4 GHz corresponds to an external magnetic 
flux in the range ±10^''<I>o which is a useful value for a 
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real experiment. However, for the RF SQUID, the flux 
required is rather less: ±3 • 10~^$o- 

One effect which is not taken into account in our sim- 
ulations is the presence of thermal photons in the cavity. 
Thermal effects can, in general, be neglected when work- 
ing with optical cavities due to the very small average 
number of photons at those frequencies. In experiment on 
solid state qubits this is, however, not generally true since 
they are operated in the microwave range. Also, the rele- 
vant temperature scale is set not only by the the phonon 
temperature (e.g. the temperature of the mixing cham- 
ber of a dilution refrigerator) but also by the amount of 
noise (essentially "hot" photons) which reaches the sys- 
tem via the leads. However, in a well-filtered system a 
total temperature of 50 mK is attainable. The resonator 
will then nearly be in its ground state with an average 
thermal occupancy h of 0.009 and thermal fluctuations in 
the photon number of the order of 0.1. This justifles ig- 
noring thermal effects in our simulations for now. That 
said, even a moderate increase in temperature can sig- 
niflcantly change the outcome of an experimental. One 
further simplifying assumption in the model is that Ti 
and T2 do not change as the qubit is detuned from the 
optimal bias point <&o/2. While this is clearly unrealistic, 
it has been experimentally shown^i that neither param- 
eter should change dramatically in the parameter range 
considered here, giving some justification to this approx- 
imation. 



V. APPLICATIONS OF CIRCUIT-QED 

One of the most important applications of circuit QED 
is the generation of single microwave photons on demand. 
Single photon sources in the optical regime have been re- 
alized using e.g. cavity QED with atoms and high finesse 
optical cavities^^. Design and fabrication of deterministic 
sources that operate in the microwave regime have proved 
to be more difficult, but a source based on superconduct- 
ing circuit-QED was recently demonstrated^^. This kind 
of source could be used for quantum radiometry, as well 
as for quantum information applications such as quan- 
tum key distribution. 

Various schemes can be envisage d^^i'^°'^"'^i'^^ , by which 
single photons can be generated with a circuit QED de- 
vice. Below, we describe a straightforward technique, 
based on manipulation of the qubit state with microwave 
pulses and rapid changes in the DC magnetic field to tune 
the qubit in and out of resonance with the cavity. 

Our technique begins with the qubit far detuned from 
the cavity and the combined system in its ground state 
|0 I). A microwave 7r-pulse is applied to the qubit to 
excite it to the state |0 t) [Fig- EU^)]; this is fol- 
lowed by a step in the magnetic field to bring the qubit 
into resonance with the cavity [Fig. EJb)]. The state of 
the system immediately after the step is still |0 t), but 
due to the qubit-cavity interaction on resonance, this is 
no longer an eigenstate, so the system begins to precess 
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FIG. 5: (Color online). Bloch sphere diagrams showing the 
state of the qubit-cavity system at successive stages of the 
single photon generation process. 



around the equator of the Bloch sphere at the vacuum 
Rabi frequency 2g. After a time 27r/4(7, the state of the 
system will be |1 |) [Fig. [51(c)]. This means that coherent 
energy exchange has taken place between the qubit and 
the cavity, creating a photon-like state. Another step in 
magnetic field detunes the qubit from the cavity so that 
the state |1 |) is once more an eigenstate of the system 
[Fig. EJd)]. The system remains in this state until the 
photon decays out of the cavity into one of the external 
waveguides in a time of order 2t:/k. By repeating this 
sequence many times, photons can be generated on de- 
mand, provided that the time window within which they 
are required is much longer than 2t:/k. 

If a scheme similar to the one above is implemented, it 
is important to prove that it generates single photons de- 
terministically, rather than stochastically. At optical fre- 
quencies this is done by studying photon-counting statis- 
tics using intcrferometric measurements. Such measure- 
ments require the use of a beamsplitter. An analogous 
experiment can be envisaged in the microwave regime, 
provided that a microwave beamsplitter can be realised. 
Such a device has been proposed recentlySS, and could 
lead to a microwave analogue of the Hanbury-Brown and 
Twiss interferometer—. 



VI. CONCLUSIONS 

We have shown that it should be possible to reach the 
strong coupling regime using a flux qubit coupled to a 
coplanar waveguide resonator. If realized, it would open 
the door to potential applications in metrology, quantum 
communication and experimental tests of quantum me- 
chanics. The fact that conventional lithography can be 
used to fabricate the samples and that the experimen- 
tal parameters can be chosen freely can in some cases 
a be significant advantage compared to CQED implen- 
tantaions utilizing e.g. atoms of ions. Since the relevant 
frequencies are in the microwave regime it is also possible 
to use well established methods to manipulate the sys- 
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tern. The main drawback compared to experiments done 
at optical frequencies is the short coherence time of the 
qubit and the fact that the system must be operated at 
very low temperatures. 
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APPENDIX A: DERIVATION OF THE 
HAMILTONIAN 



However, assuming the interaction with the environment 
is Markovian the evolution of the (reduced) density ma- 
trix of the system can be described by a master equation 
p = hp of Lindblad form^ 

k=l ^ ^ 

(Bl) 

where Cfc are Lindblad operators. In the case considered 
here we have 3 Lindblad operators. Firstly, the relaxation 
from the excited state to the ground state at a rate 71 = 
1/Ti represented by a Lindblad operator proportional 
to the lowering operator a~ , i.e. C\ = ^J^a~ . The 
cavity is loosing energy at a rate k = LOr/Q which leads 
to the "destruction" of photons in the system, C2 ~ ^/na. 
Finally, we also need to consider pure dephasing of the 
qubit at a rate 7^ = l/T^ = I/T2 - l/2Ti where T2 is 
the usual total dephasing time of the qubit. This process 

is represented by the operator C3 = \r^^z- 



It is useful to compare the informal procedure used 
in the introduction of this paper to derive the Hamil- 
tonian with a more formal approach. Starting with 
the bare qubit Hamiltonian — ^(ecr^ -I- Act^^) where e = 
2/p($a; — $0/2) we proceed just as before by noting that 
the flux threading the qubit loop will be modulated via 
the mutual inductance M that couple the fluctuations in 
the cavity to the qubit. Writing the total external flux 

as = -t-(5$, where 5$ = M^l^{a^ +a), adding 

the Hamiltonian for the oscillator mode and the external 
field and finally transforming into the eigenbasis of the 
qubit we get 

H = huir ^a^a + — ^ H — ^'^^ ^ Tig{a^(j^ + cr+a) sinO 

+^{e'''^*a^ + e*"*a) - hg{a'' + a)a^ cos 6 (Al) 

in the RWA. Here we have introduced the mixing angle 
6 = arctanA/e. This Hamiltonian is identical to the 
J-C Hamiltonian ([7]) except that we now have an effec- 
tive coupling g sin 9 and an extra term hg{a^ + a)az cos 9 
which is zero when the qubit is operated at the degener- 
acy point 9 ~ 11/2. By moving to an interaction frame 
rotating at the drive frequency oj we see that all terms in 
the Hamiltonian (|Aip are time-independent except the 
last term which picks up a factor exp(— iww), meaning 
it can be neglected in the rotating wave approxiamtion. 
Note, however, that this additional term can potentially 
play a role in the dispersive regime. 



APPENDIX C: CALCULATION OF THE 
SPECTRUM 

Our aim is to calculate the steady-state power spec- 
trum S(uj) of the intracavity field, formally this is de- 
fined in terms of the photocount output from the cavity 
as seen by a monochromatic detector-S^. The spectrum 
can be calculated from the 2-time correlation function^ 
{a\t + T)a{r)) 

1 1"°° 

S{uj)^— e"'^^a\T + t)a{t))dT (CI) 
27r J_oo 

which can be evaluated using the quantum-regression 
theorem 

< a^T + t)a{t) Tr{ate^^ap} (C2) 

where the Liovillian (which includes the three Lindblad 
operators defined in Appendix |B]) is given by the right 
hand side of equation IBll and p is the steady-state den- 
sity matrix which is the solution to Lp = These calcu- 
lations are straightforward using the built-in routines of 
the " Quantum Optics Toolbox" . 



APPENDIX B: DISSIPATION 



The effects of the environment on a quantum system is 
in general very difficult to model but is nevertheless cru- 
cial to understand since it is the cause of decoherence. 
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